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Abstract 

We discuss the problem of the stability of the isotropy of the uni- 
verse in the space of ever-expanding spatially homogeneous universes 
with a compact spatial topology. The anisotropic modes which prevent 
isotropy being asymptotically stable in Bianchi-type VII^ universes 
with non-compact topologies are excluded by topological compactness. 
Bianchi type V and type Vllh universes with compact topologies must 
be exactly isotropic. In the flat case we calculate the dynamical de- 
grees of freedom of Bianchi-type I and VIIq universes with compact 
3-spaces and show that type VIIq solutions are more general than 
type I solutions for systems with perfect fluid, although the type I 
models are more general than type VIIq in the vacuum case. For 
particular topologies the 4-velocity of any perfect fluid is required to 
be non-tilted. Various consequences for the problems of the isotropy, 
homogeneity, and flatness of the universe are discussed. 

1 Introduction 

The problem of providing a compelling explanation for the isotropy and ap- 
proximate flatness of the Universe has been the subject of extensive analysis 
ever since the discovery of the temperature isotropy of the microwave back- 
ground radiation. Historically, a variety of approaches have been taken to 
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solving these problems. The first was to justify choosing special initial con- 
ditions, at t — or t — — oo, or 'almost' initial conditions imposed at the 
Planck epoch t ~ 10~ 43 s which marks the threshold of quantum gravity. The 
second was to seek out general physical processes which might transform a 
wide range of initial conditions into a state that is similar to the presently 
observed universe after billions of years of expansion (Misner 1967). This 
strategy was pursued first within the context of the 'chaotic cosmology' pro- 
gramme during the period 1967-80. The rules usually imposed were that clas- 
sical dissipative stresses should obey the strong energy condition of general 
relativity. Typical scenarios considered involved the classical dissipation of 
anisotropy by collisional or collisionless transport processes, and the depletion 
of irregularities by the quantum particle production process (Doroshkevich, 
Zeldovich & Novikov 1968, Stewart 1968, Collins & Stewart, 1971, Zeldovich 
and Starobinskii 1972, Barrow 1977, Barrow & Matzner 1981). 

One mathematically well-defined way of approaching the question of the 
naturalness of the isotropy of the universe was to investigate the stability of 
the isotropic Friedmann universes with respects to the set of all anisotropic 
solutions of a given gravitation theory. This was first done by Collins and 
Hawking (1973) in the context of spatially homogeneous general relativis- 
ts cosmologies with zero cosmological constant by considering the stability 
of non-collapsing isotropic universes in the space of Bianchi type universes 
with the natural M 3 topology. This reduces to the study of the stability of 
particular solutions of ordinary differential equations but is made non-trivial 
by the appearance of eigenvalues with vanishing real part, and so the stabil- 
ity is generally determined at non-linear order in the expansion around the 
isotropic universes in general. 

The most general Bianchi universes containing the open Friedmann uni- 
verse as a special case are the Bianchi type Vllh spaces. The open Fried- 
mann solution was shown by Collins and Hawking not to be asymptotically 
stable in the space of all Vllh initial data if matter obeyed the strong en- 
ergy and positive density conditions. A more detailed investigation (Barrow 
1982, Barrow & Sonoda 1986) revealed that as time t — > oo isotropy is stable 
(but not asymptotically stable) and identified the attractor as a family of 
exact vacuum plane- wave spacetimes found by Lukash (1974). It contains 
the isotropic Milne model as a special case. 

The original result of Collins and Hawking is widely cited as showing that 
isotropy was unstable (see, for example, Kolb and Turner 1990). However, 
the asymptotic stability analysis of open universes is somewhat deceptive 
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because as t — > oo the imposition of the strong energy condition ensures 
that these models approach vacuum (curvature-dominated) solutions. Thus 
the asymptotic behaviour of the vacuum solutions ultimately determines the 
asymptotic stability properties of the isotropic non-vacuum solutions. How- 
ever, our past cosmological history contains very little time (if any at all) 
during which vacuum behaviour could have dominated the Hubble expan- 
sion. Thus the asymptotic stability theorems are dictated by conditions 
which have not existed in our past for any significant time interval. They 
mainly tell us about the future. Thus, even if isotropy were an asymptotically 
stable property of open sets of Bianchi Vllh initial data, it would not ex- 
plain the isotropy of our universe, because it would tell us only that isotropy 
is approached during the late curvature-dominated phase of expansion at 
redshifts less than z ~ 10. 

If the strong energy condition is violated then the conclusions will change. 
Vacuum stresses of the sort that can drive inflation must violate the strong 
energy condition in order to accelerate the mean scale factor of the universe 
and isotropy can become asymptotically stable in accordance with cosmic 
no hair theorems (Barrow 1982, Barrow & Sonoda 1986). A particular case 
of this result would be the asymptotic approach to the de Sitter metric if a 
positive cosmological constant was admitted. Violation of the strong energy 
condition is necessary but not sufficient for isotropization to occur in this 
way. 

If homogeneous initial conditions are confined to those in an open neigh- 
bourhood of the spatially flat Friedmann universe then Collins and Hawking 
(1973) went on to show that, in the space of Bianchi type VII Q initial data 
containing the flat Friedmann model as a special case, isotropy is asymptot- 
ically stable if the matter content is restricted to have zero pressure ('dust') 
to first order. The condition for isotropization requires a/H — > ast^oo, 
where a is the expansion shear and shear and H is the mean Hubble expan- 
sion rate. Isotropy is not approached during in radiation-dominated models 
of this Bianchi type. Thus if there was some reason for the initial data to 
be of zero curvature ('flat') it appeared that isotropy might be regarded as 
an asymptotically stable cosmological feature. However, again, the obser- 
vational consequences are limited because there have been so many more 
e-foldings of cosmic expansion during an era of radiation domination (from 
10~ 43 s to 10 10 s) than one dominated by dust (from 10 10 s to 10 17 s) during 
the history of the universe (Barrow 1982; Barrow & Sonoda, 1985, 1986). 
Also, as in the open universe case, if the strong energy condition and zero 
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pressure condition restrictions are relaxed to permit the inclusion of stresses 
with p + 3p < 0, where p is the matter density and p is the pressure then 
isotropy becomes asymptotically stable for an open set of initial data. This 
situation has been investigated in detail by Wainwright et al (1999) and 
Nilsson et al (2000) who discriminate between isotropization of the expan- 
sion shear, a, and that of the Weyl curvature anisotropy relative to the mean 
Hubble expansion rate, H. Radiation and dust-dominated VIIq universes 
exhibit isotropic shear isotropization, o/H — > Oast — > oo but the Weyl 
curvature anisotropy, W/H grows in the radiation dominated era and ap- 
proaches a non-zero constant (which need not be small) in the dust era as 
t — > oo (Collins and Hawking did not require the Weyl curvature anisotropy 
to tend to zero asymptotically in order for isotropization to occur in type 
VIIq). This feature was also stressed by Doroshkevich et al (1973) and 
Lukash (1974). Unlike the shear anisotropy, a large Weyl anisotropy does 
not require a large temperature anisotropy in the microwave background. 
The type VII h asymptote does have W/H — > as t — > oo because it is a 
plane-wave spacetime for which all scalar curvature invariants are zero. 

These results were interpreted as showing that if the strong energy con- 
dition holds (no inflation allowed) then generic spatially homogeneous ever- 
expanding anisotropic universes do not become isotropic during their ex- 
pansion histories. A separate analysis of closed Bianchi type IX universes 
with S 3 topology shows that isotropy is unstable under the similar condi- 
tions (Doroshkevich et al 1973). Of course, if the strong energy condition is 
dropped then a finite period of inflationary expansion in the past can reduce 
any initial anisotropy observed inside our horizon below any given level by 
a pre-specified time and drive the expansion of any open or not too closed 
universe very close to flatness for very long intervals of time at late times 
without either isotropy or flatness being approached as t — > oo. If it is too 
closed it may collapse before inflation can occur. Moreover, if further con- 
ditions of physical reality are imposed, so that the anisotropy energies in 
long-wavelength gravitational wave degrees of freedom at the Planck epoch 
do not significantly exceed the Planck density then the size of anisotropies 
today is constrained to be quite close to that observed in the microwave 
background (Barrow 1995). 

The basis of these results about the stability of isotropic universes when 
the strong energy condition holds is the instability of isotropic universes to 
perturbations by spatially homogeneous anisotropic modes. The most gen- 
eral Bianchi type universes are characterized by 4 constant parameters on 
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a Cauchy surface of constant time in vacuum and by 8 parameters in the 
presence of a perfect fluid (Ellis and MacCallum 1969). The most general 
Bianchi types containing the flat, open, and closed Friedmann universes are 
of Bianchi types VIIq, Vllh, IX, respectively and are specified by 3, 4, 4 con- 
stant parameters in vacuum and 7, 8, 8 parameters in the perfect fluid case, 
respectively For comparison, in the presence of a perfect fluid, the flat Fried- 
mann universe has no free parameters, and the open and closed models have 

1 free parameter. 

All these results assume that the topology of ever-expanding universes 
is the natural R 3 topology so that they have infinite spatial volume if the 
3-curvature is never positive. We shall show that the conclusions change 
significantly if it is assumed that the topology of the spatial sections of spa- 
tially homogeneous anisotropic universes of non-positive curvature are spa- 
tially compact. 

2 Spatially Compact Bianchi type universes 

There has been considerable interest in the possibility that the Universe 
might possess compact space sections of non-positive curvature (see the con- 
ference proceedings collection edited by Starkmanl998). These considera- 
tions have been motivated in part by the 'naturalness' of finite spatial sec- 
tions in quantum cosmologies. There is a long history of occasional inves- 
tigations of the observational consequences of compact flat universes with 
3-torus topology (Lachieze-Rey and Luminet 1995) and the possibility was 
even considered by Friedmann (1924) to show that open universes were not 
necessary spatially infinite. Studies of these compact cosmologies have been 
dominated by studies of multiple imagery (Sokolov and Shvartsman 1975) but 
there has been also some work which related the observed homogeneity and 
isotropy of the universe to its spatial topology. For example, the possibility 
of invoking a 'small' compact open or flat universe to account for the large 
scale uniformity of the universe has been proposed by Ellis and Schreiber 
(1986) following the pioneering consideration of topological effects on local 
cosmological structure by Ellis (1971). They have explored the possibility 
that the cosmological horizon problem can be significantly ameliorated by 
identifying points in such a way that observers can see all the way round 
the universe after some given time. Although this construction is easy to 
achieve in flat Friedmann universes, it is necessary to explore its effects upon 
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a fully anisotropic and inhomogeneous universe. In a related study, Ellis and 
Tavakol (1994) have considered the effects of geodesic mixing in compact 
open universes on the propagation of microwave background photons after 
last scattering and the topological constraints on anisotropy damping by this 
process were considered by Reboucas et al (1998). 

The rapid progress in observations of the cosmic microwave background 
is also providing a possibility to observe directly the topology of the universe. 
In particular, it has recently been recognized that the changed appearance of 
the microwave background sky pattern offers a sensitive probe of the spatial 
topology if the scale of periodicity is sufficiently close to the particle horizon 
scale (Levin et al, 1997, 1999, Cornish, Spergel and Starkman 1998). In the 
flat Friedmann universe there is no reason why these two scales should be 
similar. However, open universes provide a physical curvature scale which 
might be closely related to the overall periodicity scale. For this reason 
there has been much recent interest in the observational features of compact 
open Friedmann universes. Their behaviour is considerably more diverse 
than that of compact flat universes. Geodesies display chaotic divergences 
on negatively curved spaces and possible topologies are extremely complex, 
the eigenmode problem is unsolved in general, and their specifications are 
only partially understood. 

Let us now consider an extension of the problem of 'why the universe is 
isotropic?' to the situation of compact, spatially homogeneous, anisotropic 
universes. We will consider the compactification of Bianchi type universes 
with zero and negative curvature. As in the investigation made by Collins 
and Hawking (1973) for the non-compact case, we first ask what are the most 
general Bianchi types which contain the flat and open Friedmann universes as 
particular cases. This problem has been studied by several authors (Ashtekar 
and Samuel 1991, Fagundes 1985, 1992, Kodama 1998, Koike et al 1993, 
Tanimoto et al 1997, 1997a). 

If a group contains a subgroup T which acts on a manifold X as a covering 
group so that X/T becomes compact then the geometry of the manifold is said 
to admit a compact quotient. Thurston (1979) classified all maximal simply- 
connected three-dimensional geometries which admit a compact quotient into 
one of eight possible cases. 

The Bianchi classification of spatially homogeneous universes is derived 
from that of the three-dimensional Lie groups that act freely on a four- 
dimensional spacetime manifold as an isometry group with spatial orbits. 
We say that a homogeneous manifold is simply homogeneous if the group 
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of isometries has a three-dimensional subgroup that acts simply transitively 
on the manifold. The Bianchi types are subdivided into two classes (Ellis 
and MacCallum, 1969): Class A contains types / , II ,VIq,VIIq,VIII and 
IX while Class B contains types IV, V, III, Vlh and Vllh- Except for types 
IV and Vlh, we can construct a spacetime manifold with compact space for 
these models by suitable identifications of their spatially homogeneous hy- 
persurfaces. This operation will usually lower the dimension of the isometry 
group so that the spatial hypersurfaces of the resulting spacetime is no longer 
simply homogeneous. For example, it is known that a homogeneous space 
section in any Class B model cannot be simply homogeneous if it is compact. 
Compact quotients are often not homogeneous globally at all, (i.e., the full 
symmetry group does not act transitively). This happens even for the simple 
manifold T 3 j%2 which is a quotient of the Euclidean space E 3 by a discrete 
group generated by two translations in the x — y plane and a combination 
of a translation along the z-axis and a rotation by angle n around the same 
axis. Therefore, in general, the spacetime with compact space can only be 
locally homogeneous. The same argument applies to isotropy as well. 

Further new features appear in the compact cases. First, a non-trivial 
compact topology may require the geometry and matter configurations to 
have a higher symmetry than that of a simply transitive group when the 
data is lifted to the covering spacetime. Such a situation arises when the 
identification group Y cannot be contained in a simply transitive group. This 
feature has an important consequence for the isotropization problem, as we 
will see below. Second, new dynamical degrees of freedom, called 'moduli pa- 
rameters', appear. These moduli parameters describe globally non-isometric 
deformations of the geometry which preserve the local geometric structure. 
For some cosmological models, the total number of dynamical degrees of 
freedom becomes much larger than that for the non-compact case due to the 
existence of the moduli degrees of freedom. 

2.1 Open universes 

The introduction of compactness imposes a major constraint upon homo- 
geneous anisotropic universes with negative spatial curvature. The Bianchi 
types V and VI Ih contain open Friedmann universes as special cases. By in- 
specting the subgroup structure of each maximal symmetry group, one finds 
that these groups must be subgroups of the maximal symmetry correspond- 
ing to the Thurston type H 3 if they act simply transitively on constant-time 
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spatial sections £ of the covering spacetime. Hence, by Thurston's theorem, 
if the universal covering (M, g) of a locally homogeneous spacetime (M, g) 
with compact space has a symmetry group G containing a simply transitive 
subgroup of Bianchi type V or Vllh, then G must be a subgroup of the isom- 
etry group of H 3 , and so (M,g) is written as (M,g)/Y with some discrete 
group rcG. 

It is here, when determining minimal possible G, that the Class A or 
B nature of the homogeneity group becomes important. For the Class B 
Bianchi space sections, £, there exists a non- vanishing vector field v written 
as q IJ ciXj where Xi is an invariant basis with a structure constant c 1 jk, 
cj = c 1 u and qu is the component of a homogeneous spatial metric with 
respect to Xj. Since this vector has non- vanishing divergence, V-i> = g IJ ciCj, 
it cannot be invariant under the action of T. For, if it were invariant, it 
would create a well defined vector field v' with non- vanishing divergence on 
the quotient £ = E/r. But this would lead to a contradiction since 

= / d 3 x^gW ■ v' = V • v' i d 3 x^^0. 
is Jy, 

However, one finds that if G is smaller than the connected component 
of the full isometry group of H 3 , G keeps v invariant. Therefore the cover- 
ing spacetime must be invariant under a group isomorphic to the maximal 
isometry group of H 3 and be spatially homogeneous and isotropic. 

Furthermore there is no moduli freedom in the Class B case from the 
Mostow rigidity theorem (Thurston 1979, 1982), which says that two com- 
pact hyperbolic spaces are isometric up to a constant scaling if they are 
homeomorphic. Thus compact universes of Bianchi types V and Vllh ad- 
mit only an overall change of volume scaling factor. They must therefore 
be isotropic. Hence, we have the rather surprising result that compact open 
universes of type Vllh not only make isotropy an asymptotically stable prop- 
erty of the initial data set but they permit no anisotropy to be present at 
all. Recall that in the non-compact stability analyses of the open Friedmann 
models the asymptotic behaviour was approach to a particular 2-parameter 
set of anisotropic type Vllh plane-wave universes in which the ratio of the 
shear to Hubble scalar is a constant (when this ratio is zero the isotropic 
vacuum Milne universe is obtained). These anisotropic solutions are not ad- 
mitted when the space sections are compactified. Compactification requires 
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that Bianchi type initial data containing isotropic universes must be exactly 
isotropic. 

2.2 Flat universes 

The flat Friedmann universe is a special case of Bianchi types / and VIIq . In 
these Class A universes the introduction of compactification does not restrict 
the symmetry group so strongly as in the Class B cases, and anisotropy is 
permitted. However, the spatial topology still gives some weak constraints 
on the minimal symmetry of the universal covering as shown in Table [I]. 
For example, when the spatial sections are homeomorphic to T 3 , T 3 /Z 2 or 
T 3 /Z 2 x Z 2 , under an appropriate choice of coordinates, the minimal sym- 
metry of the covering data with the Bianchi I symmetry is given by R 3 xD 2 , 
while for T 3 /Z*k(k = 3,4,6) the minimal symmetry is f 3 xO(2). Here, _D 2 
is the dihedral group {1, Ri(7r), ^(tt), #3(71")} consisting of rotations by an- 
gle 71 around three orthogonal coordinate axes; 0(2) is the group generated 
by rotations around the z-axis and _D 2 , and Ax B is a semi-direct product 
of two groups A and B. Hence, the isotropization behavior of the former 
group is the same as that of the generic open type I case, but the behav- 
ior for the latter group coincides with that of axisymmetric type I model. 
However, isotropy is an asymptotically stable property for spatially compact 
Bianchi I models with perfect fluid satisfying the dominant energy condition. 
The Bianchi type I geometry does not permit the presence of non-comoving 
velocities. 

Similarly, for T 3 and T 3 /Z 2 , the minimal symmetry of the covering data 
with type VIIq symmetry coincides with the VIIq group itself, hence the 
isotropization behavior is determined by that of generic open type VIIq mod- 
els with non- vanishing generic fluid velocity. In contrast, for T 3 /Z 2 x Z 2 and 
T 3 /Zfc(fc = 3,4,6), the minimal symmetry is VII0XD2, which requires that 
the space metric be diagonal and the fluid 4-velocity thus be orthogonal to 
the spacelike hypersurfaces of constant time. Since the isotropic type VIIq 
model with radiation is already unstable against perturbations with the lat- 
ter symmetry, (Wainwright et al 1999, Nilsson et al 2000) compactification 
does not change the conclusion on the isotropization of type VIIq Bianchi 
models. 
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Symmetry Space 



Q P N m N f N N s N s (vacuum) 



R 3 xD 2 R 3 
T 3 

T 3 /Z 2 
T 3 /Z 2 x Z 2 

VII M 3 

T 3 

T 3 /Z 2 

VII xZ 2 M 3 
T 3 

T 3 /Z 2 

VII xD 2 M 3 
T 3 

T 3 /Z 2 
T 3 /Z 2 x Z 2 
T 3 /Z fc (A; = 3,4,6) 

M 3 xO(2) M 3 
T 3 

T 3 /Z 2 (aligned) 
T 3 /Z 2 (oblique) 
T 3 /Z 2 x Z 2 
T 3 /Z fc (A; = 3,4,6) 

150(3) M 3 
T 3 

T 3 /Z 2 
T 3 /Z 2 x Z 2 
T 3 /Z fc (A; = 3,4,6) 



n 
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3 
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3 





4 
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12 


10 
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10 


8 
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Q 
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11 
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1 
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1 
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8 


6 
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1 



Table 1: The number of dynamical degrees of freedom for locally homoge- 
neous systems with perfect fluid on the spaces of Thurston type E 3 . N s 
represents the dimension of the solution space. See the Appendix for the 
details of the notation. 
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Bianchi Type Vacuum Perfect fluid 

Non-compact Compact Non-compact Compact 



/ 1 10 2 11 

II 2 6 5 9 

VI 3 4 7 8 

VII 3 8 7 12 

IX - 4-8 



Table 2: Maximal degrees of freedom for spatially open and compact Bianchi 
vacuum models and for those with perfect fluid. 

2.3 Parameter Counting 

A parameter-counting classification of Bianchi type universes can be per- 
formed in the compact CclSG. clS IS explained in the Appendix in detail. Some 
of the constant parameters represent the moduli freedom, but in general 
it is difficult to distinguish clearly between the dynamical freedom and the 
moduli freedom. The general result of Kodama (1998), Koike et al (1993), 
and Tanimoto et al (1997, 1997a) is that the number of dynamical degrees 
of freedom, i.e., the dimension of the solution space, becomes larger in the 
compact cases than in the non-compact cases for the Class A Bianchi uni- 
verses. This is partly due to the appearance of the moduli degrees of freedom 
and partly due to the decrease in the freedom of diffeomorphisms connecting 
physically equivalent solutions. For example, in the case of simple 3-torus, 
T 3 , we can obtain the correct parameter count for the vacuum case by the 
following naive argument. First, in order to specify the lattice in the Eu- 
clidean 3-space so as to define the 3-torus in a rotationally invariant way, we 
need 3 parameters to specify the lengths of 3 vectors generating the lattice 
and 3 parameters to specify the relative direction angles of these vectors. 
Hence, adding their time derivatives, we need 12 parameters. If we take into 
account the Hamiltonian constraint and the time translation freedom, the 
total number reduces to 10. 

One significant feature of the parameter count in the spatially compact 
open universes is the difference between the counts for vacuum and perfect 
fluid models. If we separate the moduli and dynamical freedoms in the way 
explained in the Appendix, the diffeomorphism-invariant phase space can be, 
roughly speaking, written as the product of the phase space of the diagonal 
homogeneous system and the moduli space, at least for the vacuum Bianchi 
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/, //, VI and VIIq systems. We can show that the moduli parameters de- 
fined in this sense become constants of motion. Hence, the dynamics are 
essentially represented by the diagonal system with 3 + 3 parameters. The 
Hamiltonian constraint should then also be imposed. In this representation, 
the phase space of the non-compact case is obtained by discarding the moduli 
freedom and taking into account the additional equivalence relation for the 
diagonal system. For example, in the type / case, 3 metric coefficients of the 
diagonal system can be set to unity by a change of spatial coordinates at the 
initial time, and we obtain a final count of 6 — 3 = 3. The Hamiltonian con- 
straint and the temporal gauge fixing reduces it further to 3 — 2 = 1. In the 
compact case, we cannot do such a rescaling because it changes the moduli 
parameters, and the count is 10. From Table [I], we find that it is 8 for the 
vacuum VIIq system on T 3 . This implies that among the spatially compact 
vacuum Bianchi models, the type I model is the most generic unlike in the 
usual situation with non-compact topology, where it is the least generic. 

In contrast, when a perfect fluid is present, Bianchi type VIIq universes 
need not be diagonalizable because the momentum constraints simply relate 
the non-diagonal components of the extrinsic curvature to the spatial com- 
ponents Uj of the fluid 4-velocity. Since the momentum constraints require 
Uj to vanish for type I models, it turns out that the parameter count for 
type VIIq models is always larger than that for type I models in any given 
space topology. Therefore, the most general locally homogeneous perfect-fluid 
spacetimes that include the flat isotropic model are the Bianchi type VIIq 
models in the spatially compact case just as they are in the non-compact 
case. 

In the table we have also included, for comparison, the parameter counts 
for the compact Bianchi type IX universes with S 3 topology which con- 
tain the closed Friedmann universes as isotropic sub-cases. The vacuum IX 
universes are 4-parameter, while the perfect fluid type IX universes are 8- 
parameter. 

Here note that it is easy to extend these parameter counts to systems with 
scalar fields. Since the homogeneity-group preserving diffeomorphisms act 
trivially on scalar fields and scalar fields do not contribute to the momentum 
constraint in a locally homogeneous spacetime, the parameter count simply 
increases by two for each real component of the scalar fields irrespective of 
the Bianchi type and the other matter contents. Hence, the compact Bianchi 
I models are still more general than the compact Bianchi VII models if the 
models contain only scalar fields. 
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When counting 'parameters' in the compact cases it is important to 
remember that the lengths and the angles of the identifications are time- 
dependent variables specifying the compact spatial geometry at each time. 
Since the Einstein equation is second order in time, we must also specify their 
time derivatives, at an initial time, and then their future time develop- 
ment is completely determined. For the simple 3-torus case, six combinations 
of these variables and their time derivatives become the moduli parameters, 
which turn out to be constant in time because the Hamiltonian constraint 
does not depend on them. The other six combinations correspond to the 
diagonalized metric components and their time derivatives. However, such 
a simple counting argument sometimes fails for other topologies because the 
canonical structure becomes degenerate in the moduli sector; that is, some 
moduli parameters do not have conjugate momenta, as was shown in Kodama 
(1998). Since the existence and the uniqueness up to diffeomorphisms of the 
generic initial value problem holds for the Einstein equations, we can deter- 
mine the number of independent solutions even in such cases by calculating 
the dimension of the locally homogeneous sector of the full diffeomorphism- 
invariant phase space. This is the approach adopted in Kodama (1998), and 
differs from that adopted in the papers by Koike et al (1993), and Tanimoto 
et al, (1997, 1997a) in which spacetime solutions are classified. In the latter 
approach, knowledge on the explicit form of spacetime solutions is required. 

Finally, note that in contrast to the cases where space has R 3 topology, it 
is meaningless to compare the parameter counts for spatially compact Bianchi 
models belonging to different Thurston types such as E 3 , Nil (type II) and 
Sol ( type VIq), because spaces belonging to different Thurston types are not 
homeomorphic with each other. 

2.4 Self-similar Bianchi types 

It is of great interest to see if our results can be extended to inhomogeneous 
and anisotropic universes. In general this is a very difficult mathematical 
problem. One simple class of inhomogeneous relatives of the Bianchi type 
universes is provided by their self-similar extensions, first found by Eardley 
(1974). It can be seen that there are no self-similar spacetimes with compact 
spaces if the self-similar group is not an isometry group and acts simply 
transitively. For, as Eardley shows, in such a case the self-similarity group, 
H, contains an isometry group, G, with dimG = dimiJ — 1, and each orbit 
of G is two-dimensional. In the case of an open universe, this gives M 3 the 
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structure of a fibre bundle on a line R with fibres given by the orbits of G. If 
it can be compactified by some discrete group, T, then T must be contained in 
the isometry group G. Hence taking the quotient with Tmay compactify the 
fibres but it leaves the underlying line R unchanged. Therefore the quotient 
space is at most a fibre bundle on a line with a compact 2-space, and cannot 
be a compact 3-space. 

3 Discussion 

We have shown that the stability of open and flat isotropic universes in the 
space of spatially homogeneous initial data is strongly affected by the global 
topology of the universe. When the topology is compact, all Class B Bianchi 
type universes are forbidden unless they are isotropic. In particular, this 
means that the universes of Bianchi type Vllh, and V which contain the 
open Friedmann universes must be isotropic. 

The most general anisotropic universes containing the flat Friedmann uni- 
verses are of Bianchi type VIIq when the topology is non-compact. However, 
in the presence of a compact topology the most general flat anisotropic uni- 
verses are of Bianchi type / in the vacuum case and of type VIIq in the 
presence of a perfect fluid. Thus, in the perfect fluid case the stability prop- 
erties of the isotropic models at late times isotropic models is the same as for 
universes with non-compact topologies with non-decay of anisotropic curva- 
ture distortions as t — > oo in the radiation and dust dominated solutions, as 
found in earlier studies. 

One consequence of this new behaviour in compact spaces is to change 
the familiar conclusion that open universes are more general than flat uni- 
verses because they require more parameters for the specification of their 
initial data. When homogeneous anisotropy is present the open universes 
are forbidden while flat universes are allowed. This may have important con- 
sequences for the assessment of the significance of the 'flatness' problem in 
universes with compact topologies. It also indicates that in compact open 
and flat universes there is a close link between the properties of isotropy 
and flatness. These links should be investigated further in the context of 
inhomogeneous cosmologies. 

The local isotropy of Bianchi Vllh models with compact space implies 
that the isotropy problem is replaced by the homogeneity problem for inho- 
mogeneous perturbations of isotropic universes. That is, if there are some 
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physical processes which homogenize a perturbed Vllh compact universe 
globally, they must automatically isotropize the universe because homoge- 
neous anisotropic Vllh universes cannot exist. Alternatively, there may be 
strong restrictions on the possibility of inhomogeneous open universes with 
compact topologies. Our results also suggest the possibility that the degree 
of anisotropy of the universe is constrained by the degree of inhomogeneity if 
the universe is approximately described by a Friedmann metric with a com- 
pact space of negative curvature locally. Such constraint will apply even to 
inflationary universes and might be used to determine whether or not the 
universe is negatively curved when the present universe is almost flat due to 
inflation. 

The strong impact of compactification upon the range of possible anisotropic 
and homogeneous open universes is surprising. It shows one of the ways in 
which topological restrictions can impose significant constraints on the de- 
viation of solutions of Einstein's equations from the simplest isotropic cases 
that resemble the observed universe today. These results can be generalized 
to other gravity theories which extend general relativity by adding certain 
higher-order curvature terms to the lagrangian. These higher-order terms are 
typically negligible at late times in ever-expanding universes when the higher- 
order curvature scalars become smaller than the linear terms contributed by 
general relativity. 

One possible interpretation of our results is that the simultaneous pres- 
ence of anisotropy and spatial homogeneity is a very special combination. 
In the cases of open or flat universes with non-compact topologies, the spe- 
cialness of this type of homogeneous anisotropy is disguised by the degrees 
of freedom that are permitted for the homogeneous anisotropies. However, 
the imposition of topological compactness is extremely restrictive and makes 
some Bianchi geometries impossible in many anisotropic configurations. This 
type of restriction was also evident when rotation was introduced into closed 
compact Bianchi type IX universes (Collins and Hawking 1973a, Barrow et 
al 1985). A very strong limit on cosmic vorticity is created by attempting 
to accommodate a complex rotational dynamics in a finite positively curved 
space. The constraints on cosmic vorticity are accordingly much weaker in 
non-compact flat and open universes. 

Although the spatial compactness restricts anisotropy more weakly for 
the Bianchi models including the spatially flat Friedmann universe, some 
observable relations between anisotropy and inhomogeneity may exist for 
particular space topologies. For example, if the topology of the space is given 
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by one of T 3 /Z 2 x Z 2 and T 3 /Z k (k = 3, 4, 6), then tilting of the fluid velocity 
is forbidden if the universe is locally homogeneous. This suggests that there 
may be a constraint on the degree of tilting by the degree of inhomogeneity. 
Thus it will be interesting to investigate the relation between anisotropy and 
inhomogeneity in the framework of linear perturbation theory on spatially 
compact Bianchi models. There will also be similar restrictions on certain 
types of anisotropic stress, for example those contributed by magnetic and 
electric fields or by collisionless massless particles, in these topologies. 

It is a challenging and important problem to extend our analysis to the 
case of inhomogeneous cosmologies. We found that this cannot be done 
for the self-similar extensions of the Bianchi universes classified by Eardley 
(1974). We know that inhomogeneous solutions of Einstein's equations that 
are open or flat, with non-compact topologies, have initial data specified on a 
spacelike hypersurface of constant time by 4 arbitrary functions of 3 spatial 
coordinates in the vacuum case and by 8 arbitrary functions in the per- 
fect fluid case. In the spatially homogeneous case these arbitrary functions 
become arbitrary constants. While this suggests that the general inhomoge- 
neous solution may have parts that look locally like small perturbations of 
the spatially homogeneous models, this need not be the case for flat or open 
universes with compact topologies. Our study shows that the topology of 
the universe can impose significant restrictions upon the type of anisotropies 
that it can sustain. 
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Appendix 

In this appendix we explain how to count the number of independent param- 
eters, N s , specifying the diffeomorphism classes of solutions to the Einstein 
equations for a locally homogeneous system. 

First, we briefly review how solutions to the Einstein equations for a 
locally homogeneous system are related to those for a homogeneous system 
with a simply connected space. Let M = S x K. be a locally homogeneous 
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Table 3: The number of dynamical degrees of freedom for locally homoge- 
neous systems with perfect fluid on the spaces of Thurston type Nil. 
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Symmetry Space 
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Table 4: The number of dynamical degrees of freedom for locally homoge- 
neous systems with perfect fluid on the spaces of Thurston type Sol. J is a 
discrete transformation which permutes two of the invariant basis. 

spacetime, and M = E x R be its universal covering. The metric g and the 
matter configuration <3> on M can be lifted to M. The manifolds M and M 
are related by M = M/T, where T is a discrete group of transformations 
on M which preserve each constant time slice, E(t) = E x {t}. This lift 
(g, $) must be invariant under the action of T and, from the assumption 
of local homogeneity, T is included in a larger group, G, which acts simply 
homogeneously on M. 

In the synchronous coordinates onM = ExR9 (x, t), the action of G can 
be expressed as time-independent transformations of E. Let G' be a subgroup 
of G which acts simply transitively on E, and let x° — dt and = 1, 2, 3) 
be the invariant basis on E with respect to G'. Then the lifted data on M, 
when expressed as components with respect to the invariant basis, is specified 
by a set of functions of time X(t), which obey a set of first-order autonomous 
ordinary differential equations with four constraints, corresponding to the 
Hamiltonian constraint and the three momentum constraints. In the vacuum 
case, X(t) is given by the spatial metric components Qrj and their conjugate 
momenta P IJ . In the perfect-fluid case we must also include the fluid density 
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p and the three spatial components of the fluid velocity, uj. If G is larger than 
G', then X must have additional symmetry. Since V is included in G, these 
data are automatically invariant under T. Hence, for each solution X(t) to 
the Einstein equations on M, the pair (X(t), T) determines a solution to the 
Einstein equations on M. The functions X(t) are also uniquely determined 
by the initial data X(t Q ) at some time t = to- 
All the solutions on M are obtained in this way. However, two solutions 
on M derived from (Xi(to),Ti) and (A^^o),^) may be connected by a 
diffeomorphism on M, and so be physically equivalent. It can be shown that 
this happens when and only when the two solutions are related by a time 
translation, or the two initial data are connected as X 2 (t ) = f*X 1 (t ) and 
T 2 = fTif^ 1 by a so called homogeneity-group-preserving diffeomorphism 
(HPD) / which preserves the symmetry group G in the sense that fGf^ 1 = 
G. All the HPDs form a group, HPDG(G), which is the normalizer group of 
G in Diff(M) in mathematical terminology. 

This gauge freedom can be removed to produce a unique specification 
of each diffeomorphism class in the following way. First, we introduce the 
moduli parameters as a coordinate system of a maximal submanifold A4 in 
the space of T which is transversal to orbits of the action of HPDG(G). If 
the isotropy group H of the action of HPDG(G) is trivial, the gauge freedom 
is completely removed by this procedure. On the other hand, if the isotropy 
group is non-trivial, we further introduce a set of parameters as a coordinate 
system on a submanifold in the space of Qjj which is transversal to orbits 
of the H action. If H still has a non-trivial isotropy group, we further apply 
the same procedure to other variables, say P IJ . Eventually, we obtain a set 
of parameters specifying the diffeomorphism classes of initial data for the 
locally homogeneous system, i.e., a coordinate system for the phase space 
Fi nv (M, G). The number of parameters classifying the diffeomorphism classes 
of solutions is obtained by subtracting 1 from the dimension of Pi nv (M, G). 

This method can be also applied to the standard spatially homogeneous 
system with a simply connected space by considering the case in which F 
is trivial. For example, in the vacuum Bianchi / system, X is given by 
(Qu, P IJ ), and HPDG(IR 3 ) is given by IGL(3,M) whose element induces a 
similar transformation of the matrices Q and P. By HPDG Q can be set 
to the unit matrix. The isotropy group acts on P as 0(3). Hence P can 
be put into a diagonal form. After this reduction, the action of the residual 
HPDG on the initial data becomes trivial. Hence, by taking account of the 
Hamiltonian constraint and the time translation freedom, we find that the 
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equivalence class of the solutions are specified by a single parameter, which 
will correspond to one of the Kasner indices. 

As a procedure to determine N s for non-compact Bianchi models, our 
method is more complicated than those used by Siklos (1978) or by Ellis and 
MacCallum (1969). In their methods r inv is simply obtained as S/GL(3) 
or S/SO(3), where S is the space of components with respect to a general 
invariant basis and its structure constant, (Qu, P IJ ,C I jk), and that with 
respect to an orthonormal invariant basis, (P IJ ,C ! jk), respectively. In con- 
trast, in our method, we need information on HPDG, which is in general 
obtained by a long calculation. However, in calculating N s for compact lo- 
cally homogeneous system, this information is indispensable to determine the 
moduli degrees of freedom correctly. 

Tables |l]-|| summarize the number of dynamical degrees of freedom ob- 
tained by this procedure for locally homogeneous systems with perfect fluid 
on compact spaces with Thurston type E 3 , Nil and Sol as well as the corre- 
sponding ones for Bianchi types /, VIIq, II, and VIq on M 3 . In these tables, 
Q, P, N m and Nf denote the number of independent variables to specify the 
space metric, the extrinsic curvature, the moduli freedom and the fluid free- 
dom respectively, in the diffeomorphism invariant phase space T ; nv (M, G) 
with an invariance group of the data on a space type M. N and N s represent 
the dimensions of r inv (M, G) and of the solution space. For comparison, N s 
for vacuum systems are given on the final column where (*) implies that the 
corresponding vacuum system has a higher discrete symmetry. Note that 
the total number of dynamical degrees of freedom in this table does not take 
the Hamiltonian constraint into account, and hence is greater by 1 than the 
dimension of T mv (M,G) defined above. It should be also noted that in the 
case of Bianchi type /, the lowest symmetry is not M 3 and there are addi- 
tional discrete symmetries D 2 = {1, -Ri(vr), ii 2 (?r) , ^(t^)} corresponding to 
rotations by angle n around each coordinate axis because the fluid veloc- 
ity is obliged to be orthogonal to constant time hypersurfaces due to the 
momentum constraints. If this D 2 symmetry exists, the spacetime metric 
can be put into diagonal form with respect to the time-independent invari- 
ant basis. In contrast, if the symmetry group G does not contain D 2 , such 
time-independent diagonalization is not possible even if Q and P have three 
dynamical degrees of freedom respectively. In such cases, the diagonalization 
of the metric requires time- dependent HPDs, which produce a non-vanishing 
shift vector. 

If one allows for the time-dependent HPDs, one can always put the vari- 
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ables X into the form for the system on R 3 with the same symmetry. This 
transformation transfers some degrees of freedom in Q and P back into the 
moduli freedom and makes the moduli parameters time-dependent. This en- 
larged moduli freedom is often used as the definition of moduli freedom in the 
literature. However, the description of dynamics becomes more complicated 
using this approach. 
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